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^ : Abstract 

'T^ \ The two-body random ensemble (TBRE) for a many-body bosonic theory 

\^ ' is mapped to a problem of random polynomials on the unit interval. In this 

^ ' way one can understand the predominance of 0"*" ground states, and ana- 

Q>^ ■ lytic expressions can be derived for distributions of lowest eigenvalues, energy 

O . gaps, density of states and so forth. Recently studied nuclear spectroscopic 
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properties are addressed. 
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5 , The origins of spectroscopic properties of nuclei has received renewed attention recently in 

the context of the two-body random ensembles [0-|^. These studies provide an understanding 
of which nuclear properties are robust, depending only on the model space (1- and 2- body 
^ ' interactions), and which depend on specific strengths of interactions within the space. The 

j^ ■ 0"*" ground state is an example of a robust feature. The starting point for analyses of 

TBREs are Hamiltonians of the form H = J2k^kctck + J^ijki'^ijkmcfc'j'CkCm, where c'^,Ck 
represent boson @] or fermion [^ creation/annihilation operators for a state k, and the 
coefficients Ek, Vijkm are taken as Gaussian random variables once certain physical constraints 
are imposed, such that H commutes with the generators of total spin, isospin and so forth. 
The important distinction between the TBRE and the conventional Gaussian orthogonal 
ensemble (GOE) description of many-body Hamiltonians is that the latter does not include 
correlations between Hilbert subspaces of different quantum numbers, which are essential 
to understanding ground state and low energy spectroscopic properties of nuclei. Although 
studied for some time, very little is known analytically about the bosonic TBRE ^Q. 

We consider a bosonic model which can be treated analytically, while retaining salient 
features of more complex theories. The U{A) Vibron model |Q consists of two type of bosons, 
J'^ = O"*", 1~, and is used to describe the rotations and vibrations of diatomic molecules. In 
contrast to the fermionic U{4) problem of particles in the j = 3/2 shell, the Hilbert space 
of the bosonic theory can be arbitrarily large. We will see that this model describes many 
recently observed nuclear properties in the f/(6) TBRE 0,^. The Hamiltonian is @,0] 
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where s~^{s) and ptiP/i = —p~n) are the spherical-tensor creation (annihilation) operators 
for states with J'^ = 0+ and 1~( projection /i = 0, ±1), respectively. The square brackets 
indicate angular momentum couplings and dots scalar products. Since the matrix elements 
of the 1- and 2-body interactions are proportional to A^ and A^(A^ — 1), scaling allows 
all coefficients to be Gaussian random numbers of unit variance. The random Gaussian 
variables are grouped into a vector x = {co,C2,uo,ui,vo,es,€p). The matrix elements of H 
are well known in the vibrational basis \NnpJm), where A^ is the total number of bosons, 
Hp = 0,1,. ..,N is the number of J'" = 1^ bosons, J = np,np — 2,...,1 or is the total 
angular momentum of the many-body state ( m = — J, ..., J is omitted since it adds a trivial 
degeneracy.) We also note that even (odd) J states have only even (odd) values of Up, 
resulting in an odd/even effect depending on the choice of N. While this effect is easily 
treated, we will focus only on even A^ to simplify the presentation. 

When the dimension of the Hilbert space is large, which is typical of many-body prob- 
lems, a common approach to diagonalization is the Lanczos method. Here successive itera- 
tions of the Hamiltonian on an arbitrary trial wavefunction \E'o are performed to reduce the 
Hamiltonian to a tridiagonal form. For a given spin J of an A^ -body state, we label the 
basis hj Up = J, J + 2, ...,N — 1 or N. Clearly the dimension of the J = Hilbert space is 
maximum, while that of J = A^ is minimum. While it has been argued for some TBRE's 
that the dimensionality or the width of the lowest eigenvalue distribution is responsible for 
the preponderance of 0"*" ground states, we will see that this is not the origin here. 

We choose the trial Lanczos state to be that with np = J. Enumerating the states 
by an index k = {np — J)/2, the Hamiltonian (1) assumes the tridiagonal form H^k = 
(3k^k + Oik-i^k-i + Oik+i^k+i where ak, 13k are the matrix elements of H in this basis. 
The study of tridiagonal matrices is intimately linked to orthogonal polynomials and their 
recursion relations, where k might represent the order of the polynomial ^. In particular, the 
family of recursively generated, real polynomials are related to the characteristic polynomial 
Di[E) = det(7ij — E ■ li) of the i x i tridiagonal matrix Tij. Consequently the zeros of the 
polynomials Di are related to the eigenvalues of Tij. There are several theorems which have 
developed bounds for the zeroes of D^ as well as expressions for the extreme eigenvalues 
which we can use to derive properties of the TBRE |§|,0. 

We consider two cases here: large A^ and N = 2. For N = 2, the order of the interaction 
(2-body) is equal to the number of particles, and we expect to recover the GOE. The only 
allowed states are J = 0, 1, 2. The J = 1, 2 Hilbert spaces are 1-dimensional while J = is 
2-d. It can be readily checked that the density of J = states gives the Wigner semi-circle 
(i.e. the 2x2 matrix can be expressed as one with GOE measure), while the J = 1,2 are 
always Gaussian. The same is true for J = Jmax = N for any A^, since the Hilbert space for 
the maximum spin states is always one-dimensional so that the density of states is purely 
Gaussian (Fig. 1(b)). 

In the following, we will use the large A^ limit (typically in molecules A^ ~ 100 [^|3; 
but N = 8 calculations already agree with our analytic predictions below). We next define 
z = n„lN and 7 = J IN and construct the functions aiz, j), Biz, j) from the analytic matrix 



elements of the tridiagonal matrix, i3{z, j) = •jiz'^ + 72-2 + 73J, a{z, j) = 74-2(1 — z) + o{l/N) 
where 7^ are hnear combinations of the random coefficients x. a and j3 are the off-diagonal 
and diagonal matrix elements, respectively. Then, in the large N limit, the lowest eigenvalue 
for each spin has the form P,p!0| 

E^i„ = inf(/,(z)), f,{z)=P{z,j)-2\a{z,j)\. (2) 

In terms of the parameters of H in Eq. (1) and the matrix elements of the interactions P,^, 
we can express fj{z) = az^ + bz + dj, where z G [j, 1] and 
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b = r2 ■ X = Ep - Es + {ui - Uo - |^o|/v6) , (3) 

dj = r3- X = Es + uo/2 + f{c2 - co)/6. 

To check this approach, consider only the off-diagonal interactions {f3 = 0). Then we 
expect for small J that the minimum energy approaches £"^111 — — |'yo|A^/4v6(A^ ~ !)• The 
numerical value, denoted E^^i^, can be seen to converge to this result. For J = 0, 1,2, we 
Eiaic/ELin = 1-06, 1.05, 1.03 for N = 20; 1.01, 1.01, 1.01 for A^ = 80. Hence the desired 
spectral properties of H can be recast in terms of random polynomials on the unit interval. 
This should be generally true for any bosonic or fermionic theory since the Lanczos approach 
can be applied to either. The main effort is to determine whether there is any approximate 
analytic behavior of the functions a,P, although it is suggested to be generally true ||T0[| . 

Many general properties of random or Kac polynomials are known, but typically for 
higher order functions which are otherwise unrestricted [|ll|. To understand the properties 



of this TBRE, we first compute the distribution of coefficients. Since the variables x are 
taken as Gaussian with measure P(x) oc exp[— X]fc2;|/2], one can compute the distribution 
of coefficients of the random polynomial with V{a, b, dj) oc J dxP{x)6{a — ri ■ x)6{b — r2 ■ 
x)6{dj — r^ ■ x). Integrating yields: 

V{a, b, dj) oc exp[--J^M-^J], J = (a, b, dj). (4) 

Here M~^ = (2detM)~'^eaabef3cd{SikSji - 5ij5ki)ra,irb,j rc,krd,i is the inverse of M„^ = r^ ■ rp, 
and Ta^i is the i^^ component of r^. 

As in the U{Q) model |]1|], the frequency of a ground state 0+ is approximately 70% (for 
even N- 70.5%, 71.7%, 72.3%, 72.5% for A^ = 8,16,32,64). To understand this, we first 
compute where the minima z = Zq oi fo{z) are located. Integrating the location of the 
minima over the distribution of coefficients (4), we find the distribution of minima to be: 



^rlrl - (ri ■ r2)'^ 
2TTrjz^ + (ri ■r2)zo + rl/4' 



^(^0) = - .2j,\: :... '.2u - (5) 



In Fig. 1(a) we compare this function to results from 10^ numerical diagonalization of 
Eq. (1) for N=8,16,32 and 64 bosons. We estimate the frequency of O"*" ground states by 
first understanding where ui z G fO. 11 are the minima located. For 50% the cases, o, < 0. 



so fo{z) is an inverted parabola, and its minimum is at z = or 1, each occurring with 
equal probability (25%). For the remaining 50% of the cases, a > 0. For these cases, the 
probability of having the minimum in z G (0, 1) is 
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which has a value of 22%, leaving 14% aX z = and 14% at z = 1. Next we ask when the 0"*" 
state is a global minimum (over all J). At 2; = 0, only 0^ states are allowed (since a state 
J > has Up > J, hence z > 0). Hence for at least 39% of the cases, O"*" is the ground state. 
At z = 1, all even J are allowed, and from (2)-(3), -E^m ~ ^min = J^(c2 — co)/6. Since cq, C2 
are Gaussian random numbers, half the time 0"*" is the ground state (-E^m < ^Ininli ^-nd the 
other half, the state of maximum j, or J = N. So at z = 1, 19.5% of the ground states are 
O"*" and 19.5% are J = N. For z G (0, 1), we note that J = has even np and hence is only 
allowed at half the points. Consequently of the 22% of the minima here, no more than 11% 
can be 0"*" ground states. Hence we estimate a roughly 69.5% frequency of 0"^ ground states, 
in very good agreement with observations. 

Since it is more likely to find the minima of fj{z) on the boundary than inside, we can 
compute distribution functions of interest by restricting attention to the edges z = j and 
z = 1. While this is approximate, it does yield predictions which agree very well with the 
molecular f/(4) and nuclear ^(6) TBREs. Averaging miz{fj{z)) over P{a,b,dj) yields the 
distribution of lowest eigenvalues, denoted pj{E). These have the form 

Po{E) oc expi-E"^ /2rl)eTic{EAi) 

+Asexp{-Ey2R^)eTic{EA2) (7) 

pj_iE)^expi-Ey2R') (8) 

where R'^ = (n + ra + r^f = 1.07, Ai = (ri ■ rg + rs ■ r^)/{hr^) = 0.663, A2 = {R^ - n ■ 
r3-r2- rs)/hR = 0.632, A3 = hr^/ {Ry/ 2 det M) = 0.423, and h'^ = 2{rl{rl + 2ri ■ ra + 
^2) ~ ('"i '1^3 + ^2- f^sY), (evaluated at A^ = 64). For other values of J one can readily derive 
the general form which is a linear combination of two terms similar to (|^). In Fig. 2(a) 
we compare (^ (solid) to a Gaussian (dashes) and to results from diagonalization of H for 
selected A^. The low energy excess is readily described, and the results are not Gaussian for 
any A^, in contrast to results from the dilute limit [Q. In 2(b) we compare the same function 
with slightly modified parameters (rs = a/2 and Ai = 3/4) and see that this functional 
form can readily account for the observed asymmetry in the U{6) results [|I[]. The results 
for maximum J are shown in Fig. 1(b). Since this Hilbert space is 1-d, the distribution also 
corresponds to the level density. 

It is interesting to use fj{z) to estimate the level density pj{E) for states of spin J by 
averaging this over the space of random polynomials. For J = Jmax, we trivially recover (8). 
For J = we find 
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where g{z) = Z^MZ, Z = (z^,z, 1), which agrees well with calculations (Fig. 3). Po{E) 
has moments: 



on fl 

{E^^)j=o = -r{n + 1/2) / g{zrdz. (10) 
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One can see that the shape of the level density is a superposition of Gaussians of varying 
width. The width of the O"*" density of states is {E'^)j=q = 0.784, while for the maximum 
spin {E'^)j=N = 1-137. In this model we see that there is no direct relation between the 
widths of the distributions and the probability of having a 0"^ ground state, as has been 
conjectured in other TBREs 0,|l[. Certainly the functional dependence of these moments 
are distinct from the dependence of (5)- (6), which is related to the ground state problem. 
Hence, the level densities do not reflect ground state properties in the sense that correlations 
between subspaces of different J, which are central to that question, are not reflected in these 
functions. 

The distribution of 1^ — 0^ energy gaps, g = -E^jn ~ ^min^ denoted p{g), is obtained by 
averaging g over Eq. (4): 

M oc exp{-gy4rl) [b^ exp{-g\rl + n ■ r2)/A) (11) 

xeTic{EB2) - eiicigBs)] 

where we flnd the scaling g = ^/2Ng. Here, Bi = r2/Jr2 + 4ri ■ r2 + 4rf = 1.059,i?2 = 

(2rf + ri ■ ra - A)/2VdA = -0.03, B3 = {n ■ r^ + r|)/(2r2V^) = -0.311, d = r\rl - (n ■ ra)^ 
and A = 4(r^ + r\ ■ r-i) + r\ (the numerical value is for A^ = 64). This scaling in A^ is evident 
in Fig. 4(a) when we plot versus g. We note that p((7) describes the shape away from the 
origin. (Near the origin there is an abundance of small gaps that arise from the omitted 
region z G (0, 1)). We use the same function in Fig. 4(b) to compare to the E^+ — Eq+ gaps 
computed in the IBM, with g = Ng. The same scaling is apparent. Finally, the distribution 
of R4/2 = (£'4+ — Eq+)/{E2+ — Eq+) has been measured in the ^(6) TBRE. The analogous 
quantity in this model is R2/1 = {E2+ — Eq+)/{Ei- — Eq+). The distribution P(i?2/i) can 
be derived and has peaks at 2 and 3 (= (J + 1)/(J — 1) for J = 2), consistent with those 
in the f/(6) model. To leading order this function is just a sum of two delta functions, and 
higher order corrections must be included to get the shape. 

We have examined a bosonic TBRE which shares the salient features of more complex 
TBREs, such as a Wigner limit for small A^, Gaussian level densities for certain states, 
preponderance of O"*" ground states, and so forth. By mapping the TBRE onto random 
polynomials on the unit interval, we are able to analytically understand many properties of 
the TBRE found numerically, including the frequency of O"*" ground states. We flnd the latter 
is not attributed to the width of the level densities or the dimension of the Hilbert spaces. 
Rather, the various interactions in H tend to put the extreme values of spin (J = 0, A^) 
at the ends of the spectra, enhancing their chances to be the ground state. These results 
provide the flrst analytic understanding of ground state properties, distributions of lowest 
eigenvalues, gaps and level densities for the TBRE. These functions are also found to describe 
the nuclear properties obtained in the IBM. Scaling behavior has also been predicted and 
verifled in certain observables. Since the analytic Lanczos approach can be generally applied, 
it would be interesting to see if a more general connection can be made between the TBREs 
and random polynomials |T^ . 

I thank A. Frank, R. Bijker, R. Casten and F. lachello for many interesting discussions. 
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FIG. 1. (a) Theoretical distribution (5) of extrema of /o(^) (solid) compared to numerical 
TBRE results (symbols), (b) Theoretical eigenvalue distribution for highest spin states pn{E) 
(solid; Eq. (8)) are seen to be Gaussian and agree with numerical results (symbols). This is also 
the level density pn{E) oi J = N states for any N . (In this and the following figures, the vertical 
scale is arbitrary). 




FIG. 2. (a) Lowest eigenvalue distribution (7) for J'^ = 0+ states (solid) compared to numerical 
TBRE calculations. The Gaussian is for reference, (b) Same as (a) but for the nuclear case. 
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. Theoretical density of states (9) for 0"'" states (solid) compared to a Gaussian, and 
TBRE calculations for selected N (symbols). 
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FIG. 4. (a) Distribution of energy gaps p((/) where g = V^N{Emin{J = 1) — Emin{J = 0)). 
The description is good, demonstrating a predicted scaling, with the exception of the behavior at 
the origin, (b) Same as (a) but for the nuclear case. Here g = N[Emin{J = 2) — Emin{J = 0)). 



